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Abstract —In this paper, we investigate optimal resonrce alloca¬ 
tion in a power beacon-assisted wireless-powered communication 
network (PB-WPCN), which consists of a set of hybrid access 
point (AP)-source pairs and a power beacon (PB). Each source, 
which has no embedded power supply, first harvests energy 
from its associated AP and/or the PB in the downlink (DL) and 
then uses the harvested energy to transmit information to its 
AP in the uplink (UL). We consider both cooperative and non- 
cooperative scenarios based on whether the PB is cooperative 
with the APs or not. For the cooperative scenario, we formulate 
a social welfare maximization problem to maximize the weighted 
sum-throughput of all AP-source pairs, which is subsequently 
solved by a water-filling based distributed algorithm. In the non- 
cooperative scenario, all the APs and the PB are assumed to be 
rational and self-interested such that incentives from each AP 
are needed for the PB to provide wireless charging service. We 
then formulate an auction game and propose an auction based 
distributed algorithm by considering the PB as the auctioneer and 
the APs as the bidders. Finally, numerical results are performed 
to validate the convergence of both the proposed algorithms and 
demonstrate the impacts of various system parameters. 

Index Terms —Power beacon-assisted wireless-powered com¬ 
munication network (PB-WPCN), wireless energy transfer 
(WET), resource allocation, auction theory. 

I. Introduction 

Recently, wireless energy transfer (WET) technologies have 
drawn wide attention with their capability of energy supply Q. 
Conventionally, the very limited energy of wireless devices 
powered by batteries largely constrains the communication 
performance in many practical cases such as wireless sensor 
networks Q. Also, the battery replacement for wireless de¬ 
vices is not always convenient or feasible in many applications. 
These challenging issues have boosted the development of 
WET technologies, which enable wireless devices to harvest 
energy from wireless signals for their intended applications. As 
one category of existing WET techniques, the radio frequency 
(RE)-enabled WET 0 considered in this paper, provides the 
feasibility of a long-range energy transfer (up to tens of 
meters 0) compared to other technologies, such as inductive 
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coupling 1 ^ and magnetic resonance coupling 0 - RE-enabled 
WET has not been widely used in practice largely due to the 
high propagation loss of RE signals. However, due to the lat¬ 
est breakthroughs in wireless communications, namely small 
cells 0, transmission using large-scale antenna arrays (i.e., 
massive MIMO) 0, millimeter-wave communications tn. 
the transmission distances would be dramatically reduced and 
sharp beamforming could be enabled, which will significantly 
reduce the propagation loss and achieve much higher WET 
efficiencies 0 - Eurthermore, the energy consumption of users 
will be continuously reduced by the advancements in low- 
power electronics 0- Thus, we believe that the RE-enabled 
WET has a great potential to be widely implemented in the 
next-generation wireless communication systems. Then, by 
considering the RE-enabled WET technique, a fully wireless- 
powered communication network (WPCN) can be established 
with no need of battery replacement 0. In a WPCN, wireless 
devices are only powered by WET in the downlink (DL) and 
transmit their information using the harvested energy in the 
uplink (UL) GD- 

In open literature, there have been several papers that 
focused on the design of WPCNs for different setups 0- 
p3| . In 110 , a “harvest-then-transmit” protocol was proposed 
for a multi-user WPCN, where users first harvest energy from 
RE signals broadcast by a single antenna hybrid access point 
(AP) in the DL and then transmit information to the AP in the 
UL via time division multiple access (TDMA). Moreover, the 
DL WET time from the AP and UL information transmission 
time of individual users were jointly optimized to maximize 
the system sum-throughput. 0 extended 0 to a multi¬ 
antenna WPCN scenario, where a multi-antenna AP enables 
the simultaneous UL transmission via space division multiple 
access (SDMA). (14| and 0 considered full-duplex WPCNs, 
where a full-duplex AP is adopted to provide the simultaneous 
DL WET and UL information transmission. Moreover, 0- 
0 focused on the development of cooperative protocols for 
WPCNs with different setups. In all aforementioned papers, 
only AP is considered as the energy source of the whole 
network. In 0, the authors proposed the idea of deploying 
dedicated power nodes, named power beacons (PBs), to enable 
WET in the DL. By resorting to the stochastic geometry 
theory, the densities and transmit power of PBs are investigated 
under data links’ outage constraint. With this PB-based WET, 
we thus could consider a new network setup, named “PB- 
assisted WPCN (PB-WPCN)”, in which each user can harvest 
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wireless energy not only from the AP but also from the 
deployed PB. For this new model, a natural question that arises 
is “how to optimally allocate the resources of PB-WPCNs, 
including the PBs’ energy, and the time for DL WET as well 
as the UL information transmission?”. To the best of our 
knowledge, this is still an open question, which motivates 
this paper. 

In this paper, we consider a PB-WPCN consisting of one 
multi-antenna PB and multiple single-antenna AP-source pairs. 
In view of the state of art and trend of RF energy transfer 
0^ the considered network setup is very likely to hnd its 
applications in the practical scenario of small cells, such as 
picocells (range from 10 to 100 meters) and femtocells (WiFi- 
like range), which has been regarded as one of the key enabling 
technologies of the upcoming 5G cellular networks Q. It 
is also worth mentioning that although introducing the PB 
may result in some extra cost and complexity to the system, 
this could be benehcial as a whole based on the following 
considerations: (1) The PB could be dedicated designed for 
WET only and thus can achieve a higher WET efficiency 
by exploiting the benehts of energy beamforming enabled 
by multiple antennas |T9), j^. (2) The deployment of the 
PB could be more flexible since it has much looser backhaul 
requirement. (3) By using both the AP and PB to perform 
WET simultaneously, the transmit diversity could be achieved 
and thus could make the RE energy harvesting at the user 
side more robust. We consider that the APs and the PB are 
connected to constant power supplies. Each AP aims to collect 
the information from its associated source. It is assumed that 
each source has no embedded energy supply but has the ability 
to harvest and store the energy from RE signals broadcast by 
its AP and/or the PB. The PB is installed to assist APs during 
the DL WET phase by providing wireless charging service. 
Note that the energy allocation of the PB and time allocation 
of each AP-source pair are tangled together. This makes the 
distributed and optimal resource allocation for the considered 
PB-WPCN non-trivial at all. The main contributions of this 
paper are two-fold. 

Firstly, we consider a cooperative scenario, where all the 
AP-source pairs and the PB are cooperative to maximize the 
network social welfare, which is defined as the weighted sum- 
throughput of all AP-source pairs. This is done by jointly 
optimizing the DL WET time of each AP and the energy 
allocation of the PB. A water-filling based algorithm is sub¬ 
sequently proposed to optimally solve the formulated social 
welfare maximization problem in a distributed manner. Note 
that this cooperative scenario can correspond to the situation 
that APs and the PB are deployed by the same operator, as 
the one considered in El. 

Secondly, we consider a non-cooperative scenario, where 
all the APs and the PB are assumed to be rational and self- 
interested. This scenario could be used to model the case that 
the APs and the PB are installed by different operators. The 
PB will ask a monetary payment on its wireless charging 
service and the APs will measure its benefits obtained from 
the achievable throughput with its payment to the PB. The 
conflicting interactions among the APs and the PB can be 
modeled by game theory. Specifically, we formulate an auction 



Fig. 1. System model for the proposed PB-WPCN. 


game for this non-cooperative scenario based on the well- 
known Ausubel auction which induces truthful bidding 
and achieves the maximum social welfare |[2T1-|[24|. In this 
formulated game, the APs are bidders and the PB is the 
seller as well as the auctioneer. An auction based distributed 
algorithm is proposed to analyze the formulated game and 
its convergence is subsequently proved. Note that auction 
theory p2) has been well investigated and widely applied 
for designing the resource allocation in cognitive radio net¬ 
works p5) , D2D communication networks | |26) , cooperative 
communication networks as well as energy harvesting 
networks | |28) , p9) . A cooperative network with multiple 
source-destination pairs and a relay was considered in | [^ , 
where an auction based power allocation scheme was proposed 
to allocate the harvested energy of the relay in a distributed 
manner. 1291 formulated a non-cooperative game to model the 
competitive WET bidding of the users in a WPCN with one AP 
and multiple users, where the AP adopted auction mechanism 
for DL WET. Both |28| and adopted the concept of 


Nash equilibrium to evaluate the strategic interactions between 
bidders only. In contrast, in this paper we adopt a different 
auction mechanism to study the hierarchal interaction between 
the PB (auctioneer) and the AP-source pairs (bidders). 


The rest of this paper is organized as follows. We first 
elaborate the system model in Section II. Then, a social 
welfare maximization problem for the cooperative scenario 
is formulated and a water-filling based distributed algorithm 
is proposed in Section III. In section IV, an auction game 
for the non-cooperative scenario is designed and an auction 
based distributed algorithm is proposed. Section V presents 
the numerical results. Einally, Section VI concludes the paper. 


IT System Model 

In this paper, as shown in Pig. we consider a PB-WPCN 
with one PB and N APs associated with N (information) 
sources. Each AP aims to collect the information from its 
associated source. We denote by JV = {I,-- - ,N} the set 
of AP-source pairs and the ith AP-source pair consists of the 
ith AP and the ith (associated) source, Vi G Af. We consider 
that each source has no fixed energy supplies and thus needs 
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to replenish energy from the signa(^ sent by its AP and/or 
the PB. The APs and the PB are connected to constant power 
supplies. We also consider that all APs are connected with 
the PB via backhaul but they are not connected with each 
other directly. We assume that all the AP-source pairs work 
on orthogonal frequency bands, while the PB can work on all 
frequency bands. All AP-source pairs work in the half-duplex 
mode. The APs and sources are each equipped with one single 
antenna and the PB is equipped with M > 1 antennas. Full 
channel state information (CSI) is assumed to be available at 
the transmitter side. 

We exploit the “harvest-then-transmit” protocol proposed 
in GD- Specifically, each source first harvests energy from 
RF signals broadcast by its associated AP and/or the PB in 
the DL and then uses the harvested energy for information 
transmission in the UL. By considering that the sources may be 
low-cost, low-complexity and low-energy devices, we assume 
that an integrated architecture is adopted at each source, 
where the energy harvesting component and the information 
processing component are integrated together by using one 
rectifier circuit m- In this context, each source can only 
harvest energy from the RF signals with the frequency inside 
its working bancj^ which means that each source can only 
harvest energy from its own AP and/or the PB. 

As illustrated in Fig. 2, we assume that all APs and PB 
perform WET to sources simultaneously from the beginning 
of each transmission block T. We denote by the fraction of 
a transmission block T for DL WET from the ith AP to the 
ith source, and we denote by the fraction for DL WET from 
the PB to the ith source. Then the DL WET time for the ith 
source is thus given by max (r^, r') T- During the remaining 
time (1 — max (ri, r')) T, the ith source uses the harvested 
energy to transmit its information to its corresponding AP in 
the UL. Eor convenience and without loss of generality, we 
assume that T = 1 in the sequel of this paper. 

We assume that the DL and UL channels are reciprocal. 
Let scalar denote the complex channel between the zth AP 
and the ith source and let S (jgjjote the complex 

channel vector between the PB and the ith source, where 
^Mxi (jenotes a set of all complex vectors of size M x 1. 
Besides, we use S ([]Mxi jq denote the beamforming 
vector at the PB applied to the energy signal transmitted to 
the ith source with |jvvi|p = 1, where |jwi|| is the Euclidean 
norm of w^. It is assumed that the ith AP transmits with 
power Pi to the ith source and the PB transmits to each 

'Note that the energy signal could be designed as a zero-mean pseudo¬ 
random signal with arbitrary distribution as long as its power spectral density 
satisfies certain regulations on radio signal radiation for the operating band 
of interest since it does not carry any intentional information [30] . 

^Note that if the sources are advanced enough such that they ai'e equipped 
with separated front-end hai'dwai'e for the energy harvesting unit and the 
information processing unit, it is also interesting to consider the case that 
the RF energy transfer is performed in a wideband manner and the sources 
can haiwest energy over all frequency bands. However, the design of such a 
system with wideband energy transfer would be quite different from and much 
more complicated than the considered one with narrowband energy transfer. 
Thus, as the initial effort towards the design of the PB-WPCN, we choose 
to focus on the case with narrowband energy transfer in this paper and we 
would like to consider the design of another interesting setup with wideband 
energy transfer as our future work. 
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Fig. 2. An illustration of the time diagram for ith AP-source pair in a 
transmission block T. 


source on different frequency band^with the same power pi,. 
Eurthermore, the PB is constrained by a total energy during 
each transmission block, denoted by which leads to a 

total energy constrainj^^.g^ r/pt < El°*. 

Then, during the DL WET phase, the received signal at the 
ith source, denoted by pi, can be expressed as follows when 
both the ith AP and the PB perform WET to the ith source, 

Vi = y/Pi9iXi,AP + y/Pb^f^iXi^PB + tfi, (1) 


where Xi^AP and Xi^ps are the transmitted signals from 
the ith AP and the PB to the ith source, respectively, with 

\xi pb\^ = 1- E [•] denotes the 


E 


\XiAp\ 


= 1 and E 


expectation. (•)*^ is the Flermitian transpose, is the additive 
white Gaussian noise (AWGN) with a zero mean and variance 
(T^. Note that the first or the second term on the right-hand 
side of ([T]| should be removed when only the PB or the ith 
AP transfers energy to the ith source. 

Since the multi-antenna PB transmits energy signal to a 
single-antenna source and full CSI is available at the PB, 
the optimal energy beamforming vector at the PB should be 
maximum ratio transmission (MRT) | |20l , 1^ . We thus have 
Wi = u^- Therefore, at the end of WET phase, the amount 
of energy harvested by the ith source, denoted by Ef, can be 
written by 


Et = rj {npiGi + TiPbKi ), ( 2 ) 

where 0 < p < 1 is the energy conversion efficiency, Gi = 
I Pip is the channel power gain between the ith AP and the 
ith source, Ki = ||ki|p is the equivalent channel power gain 
between the PB and the ith source, and the harvested energy 
from the noise is ignored since it is negligible in practice 
i jiTI . We also define that G = [Gi,--- ,GAr]^ and K = 
[Ki , • • • , respectively. Note that the derivation from 

^Multi-input single-output and orthogonal frequency division multiplexing 
(MISO-OFDM) can be adopted at the PB to power the sources on different 
frequency bands via multiple antennas. 

Generally, an energy consumption constraint should also be imposed to 
each AP as that to the PB. But, the energy consumption requirement at the 
APs could be easily satisfied because we consider that the APs are connected 
to constant power supplies and only serve their own source with fixed transmit 
power. In contrast, the PB needs to transmit wireless energy to multiple 
sources and its energy consumption constraint could be frequently violated 
when the number of sources becomes large. Moreover, it is worth mentioning 
that the proposed distributed algorithm elaborated in Section III can readily 
be extended to solve the problem with additional energy constraints for the 
APs. For the purpose of exposition, in this paper we choose to ignore the 
energy constraint at each AP. 
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to Q is carried out based on the assumption that the signals 
transmitted by the AP and PB are mutually independent as in 



After the source replenishes its energy in the DL, it then 
transmits its information to the AP in the UL. It is assumed 
that the harvested energy is exhausted by the source for 
information transmission. The effect of energy storage and 
energy consumption of the circuit is disregarded for clarity as 
in Og, |Tg. The transmit power of the ith source, denoted 
by Qi, is thus equal to 


Ef ^ 7? {TiP^Gj + TlpbKj) 

1 - max{Ti, r/} 1 - max (n, r') 


(3) 


Then, the achievable throughput at the ith AP can be written as 


R^ (a, t') = (1 - max (n, t')) W loga 
= (1 - max(ri,r')) W 



X log2 


A G^r]{Tip^G^ +TlpbK,) \ 

\ (1 - max (Ti,r/)) 0-2 J ’ 


(4) 


where W is the bandwidth and is the noise power, which 
are assumed to be the same for all AP-source pairs, without 
loss of generality. It can be observed from 0 that the 
achievable throughput for the ith source to its AP can be 
increased with the assistance of the PB. 

In this paper, we consider both cooperative and non- 
cooperative scenarios based on whether the PB is cooperative 
with the APs or not. In the cooperative scenario, the PB and 
the APs cooperate to maximize the network social welfare, 
dehned as the weighted sum-throughput of all AP-source pairs. 
In contrast, the PB and the APs are considered to be self- 
interested in the non-cooperative scenario. More specihcally, 
the PB will ask a monetary compensation for its wireless 
charging service and each AP-source pair will value its benehts 
with its payment to the PBs. We will design the resource 
allocation schemes for these two scenarios in the following 
sections, respectively. 


III. Cooperative Scenario 

In this section, we hrst formulate a social welfare maxi¬ 
mization problem for the cooperative scenario and then solve 
it with a water-hlling based algorithm in a distributed manner. 


A. Problem Formulation 

We consider a weight > 0 for the ith AP, which 
represents a gain per unit throughput from the ith source to 
its AP. Then, the social welfare maximization problem can be 
formulated as 

max XiRi{Ti,Ti), 

s-t. Ti,Tl > 0, 0 < max (ri,T') < 1, Vi S Af, (5) 

iGTV 

where the constraint 0 < max(Ti,r') < 1 guarantees that the 
weighted sum-throughput is larger than zero. By observing 


the objective function of problem Q, we obtain the following 
lemma; 


Lemma 1 . The function Ri (r^, t') in 0 can be rewritten as 


Ri (a , Ei ) 


(l-r,)ITlog2 


/ Gip {npiGt + E^Ki) \ 

\ (1-A)a2 j’ 

( 6 ) 


where Ei = r^pt is the amount of energy that the PB allocates 
to the ith source. 


Proof: We hrst prove that max{ri,T'} = by contradic¬ 
tion. If A < t[, then max{Ti,r'} = r/ and 0 becomes 




(1 - r;)kPlog2 


Gip {npiGi -f T'iPbKi) \ 

(l-^)a2^‘" J- 
(7) 


We can easily see that the right-hand side of 0 is a monoton- 
ically increasing function of r^. This means that the achievable 
throughput Ri can always be enhanced by increasing the value 
of Ti to that of rf This contradicts the assumption that n < t[. 
We thus have Ti > i.e., max{Ti, r'} = Ti. Furthermore, for 

the ease of presentation, we dehne that Ei = T^pb, which is 
the amount of energy that the PB allocates to the ith source. 
Applying these two operations into 0, we can rewrite it as 
0, which completes this proof. ■ 

Based on Lemma we can reformulate the problem 0 as 

max XiRi {n, Ei), 
i&N 

s.t. 0 ^ X ^ 1, 

( 8 ) 

0 f E f rpb, 

Y,E,< El^\ 

ieM 

where t = [ri,--- ,tn]^ and E = [Ei,--- ^En]^, the 
symbols ^ and ^ represent the element-wise inequality, 0 or 
1 is a vector of zeros or ones that has the same size as t and 
E, and the constraint of E is derived based on Ei = T[pb and 
0 < r' < Ti. Moreover, we denote by {t*,E*) the optimal 
solution to the problem 0. 

We can see from 0 that the social welfare can only be 
maximized by jointly designing the DL WET time of each 
AP and the energy allocation of the PB. It is worth noting 
that T and E are mutually interdependent and the achievable 
throughput of all APs are coupled together due to the total en¬ 
ergy constraint J2ieAf Ei < El°\ 


B. Optimal Solution and Distributed Algorithm 

The convexity of the objective function in problem 0 
can readily be checked by introducing a new optimization 
variable 9i = \ — max{ri,r'}, Vi G N. Substituting 6i into 
0, we can observe that Ri {Ti,Tl) is jointly concave on Ti, 
t[ and 9 i, since 9 i is the perspective variable of the log. 
Accordingly, the problem 0 is convex and could be solved by 
applying standard convex optimization approaches. However, 
these approaches are normally done in a centralized manner. 
In practice, a distributed approach is of more interest because 
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it can significantly reduce the network overhead. Motivated 
by this, we propose a distributed method with three steps to 
resolve it in this section. In step 1, we first hnd the optimal 
relationship between and Ei by expressing Ti as a function 
of Ei. Then, the problem (j^ can be reformulated as a problem 
with only E as a variable. In step 2, we investigate the 
properties of the reformulated problem. Finally, a water-hlling 
based algorithm is proposed to find the optimal solution of the 
new problem in step 3. We describe the details of these three 
steps in the following subsections. 

1) Step 1. Problem Reformulation: We first derive the 
expression of Ti as a function of Ei for each AP. With a 
given E that satishes 0 f E ^ Iph and 
the problem ([8]l is decoupled and we can have the following 
optimization problem for each AP regarding Ti, 

Ei 

max Si {n) , s.t. 0 <Ti <1, Ti> (9) 

n pb 

where Si (Ti) can be expressed by 

(a) — ^iRi (fzjF/j) 

= XtW (1 - Ti) log 2 ( 1 + 


G,r] (Tip^Gj + EjKi) 
(1 - Ti)a‘^ 


We denote by Ti{Ei) the optimal solution of (|^, which i 
given in the following proposition. 


Hence, we now can reformulate the problem ([^ with only 
E as the variable as follows 


max 

E 


'y Si {Ei) , s.t. 0 f E ^ ^ ^ Ei < E^ 


tot 

b ■ 


ieAf 




(15) 


Note that we can solve the problem (15i to obtain E , and 
then we can calculate x* based on ( 111 , thereby solving the 
original problem ([^. 

2 ) Step 2. Property Characterization: To solve the problem 
we hrst characterize the properties of the objective 
function J2iej\r •^UEi) is a positive 

summation of N independent functions with the same struc¬ 
ture, we then first investigate the property of any function 
Si {Ei), which is summarized in the following proposition: 

Proposition 2. When 0 < Ei < pb, Si {Ef is differentiable 
and the gradient of Si {Ef, denoted by SI Si {Ef), is continu¬ 
ous. VSi {Ei) is derived as 


VS, {Ei) = 


if 0<E,< Ef^, 

< Ei < Pb, 


( 10 ) 
is 


Proposition 1. Given Ei € [0,pb), the optimal solution Ti {Ei) 
to the problem ^ can be expressed by 

( { 4 -l)a--G,vE,K, n < 

T, {Ei) = \ ^ ^ 

if, if Ef”^ <E,<pb, 

where 

Pb (4 - l) 0-2 

F;'™ = ^^-, (12) 


^4^{Ei), if El 
where oci is a constant given by 

XiWG,pK, 

CXj, — j. 5 

z] (T^ In 2 

and j3i {Ei) is strictly decreasing, which is expressed as 

X,E, 


(16) 


(17) 


Pi {Ei) = - log 2 


Pb 


1 + 


Pb - Ei 


X,WX, 


(18) 


(4 - l) cr2 -p Gip {piGi+pbKf) 
and zp > 1 can be expressed by 

4 =(w { °sZai) ) + 0 • 

in which W {x) is the Lambert W function that is the solution 
to the equality x = W exp(W). 

Proof: See Appendix A. ■ 

Then, by replacing Ti in Si (t^) with {Ef) given in 
Proposition we have a function Si {Ef with only Ei as 
the variable given by the following Lemma. 

Lemma 2. Si {Ef) can be expressed by 

S, {Ef = 

if 0<E,< Ef^, 

(l - f) (l + i^) > ‘f < Pb, 

(14) 

where X, 4 _ 

Proof: See Appendix B. ■ 


{Pb — Ei -\- XiEf In 2 

Proof: See Appendix B. ■ 

Therefore, based on Proposition 2, it is easy to show that 
Si {Ef is a concave function of Ei when 0 < Ei < pb- Then, 
by letting the first-order derivative of Si {Ef be equal to zero, 
i.e.. Pi {Ef — 0, we obtain the stationary point that maximizes 
Si {Ef, denoted by E°, given by 


E° = 


Pb - l) cr 


^z| - 1^ 0-2 _|_ {piGi -PpbKf 


where zj > 1 can be expressed as 

t_4 a, f {PrGi + PbKf - (T^ 


-hi 


(19) 


( 20 ) 


=■=“!■ Tv 

Note that E° has the same structure as Ef'^ given in (12 1 . The 
only differences between them are the parameters zl and zj. 
Furthermore, we can easily observe that zf > zj > 1 based 
on the property of the equation ( |4T] i in Appendix A. Let us 
divide both the top and the bottom of the right-hand side in 
(191 by — I'j ( in (12 1 by a^). We thus can 

obtain that Ef"^ < E° < Pb- Now, we have an important 
property of Si {Ef: it is increasing when 0 < Ei < Ef and 
decreasing when Ef < Ei < Pb- The PB will at most allocate 
the amount of energy that is equal to Ef to the ith AP to 
achieve the maximization of the social welfare. We thus only 






















6 


need to focus on the interval 0 < Ei < E° for each AR We 
now can update the property of the function Si {Ei) over the 
interval 0 < Ei < E° in the following proposition: 

Proposition 3. For 0 < Ei < E°, Si {Ei) is differentiable, 
increasing and concave, andVSi {Ei) is continuous. VSi {Ei) 
is constant when 0 < Ei < while strictly decreasing 

when < Ei < E° and equal to zero when Ei = E°. 

For the ease of understanding, functions {Ei), Si {Ei) 
and \7Si {Ei), given in ( 111 , ( [l4] i and ( |T6| ), respectively, are 
graphically interpreted in Fig. We use the black and red 
curves to characterize the functions in the intervals 0 < Ei < 


and 

that 


E:, 

E: 

The harvesting time from the PB 


Pb 


curve in the interval E: 


lim 


< E, 


max 

E 


s.t. 0 ^ E ^ E°, 

Y.e,<eI°\ 


0 has a trivial solution that E* = E° 
only consider the case that El°* < 




Et = 


^tot 

A 


-E 


jGW\{i} Ej ; 


if 0 < 1/ < cti, 

if V = tti 

if V > tti, 


where z\ > 1 is the unique solution of the equation 


i In {zi 


/ i/pb In 2 

V A,IF 


- 1 U* + 1 = 


G^r] {p^Gj +pbKi) 


(24) 


^ < Ei < Pb, respectively. We can see from 
Ti {Ei) linearly decreases in the interval 0 < 
Ef"^, but linearly increases when Eli *"*<£;*< Pb. 

t' = ^ is also shown in 


Fig. 3(a) using the blue dashed line, which is always linearly 
increasing in the whole interval and coincides with the red 


< Pb. It can be observed 
in Fig. 3(b) that Si {Ei) is linearly increasing in the interval 
0 < Ej < which corresponds to the observation in 

Fig. 3(c) that VSi {Ei) is constant for Q < Ei < E^. With 


the increasing of Ei, Si {Ei) is maximized at the point E°, 
where {Ei) equals zero. Finally, when E° < Ei < pb. 
Si {Ei) keeps decreasing and {Ei) is shown to be always 
negative. 

3) Step 3. A Distributed Algorithm: Based on the analysis 
in step 2, the optimization problem ([T5]l can be written as 


( 21 ) 


where E° = [E° • • • , Eff] . Based on the property of Si {Ei) 
summarized in Proposition we then can solve the problem 
( |2 T] i using a water-filling based approach p5| . It should be 
noted that when Eff^ > Y^i^M Ef, the optimization problem 

Vi G Af. We thus 
E°, namely, the 


Proof: See Appendix C. ■ 

According to Proposition and Proposition there exists 
a unique v to meet the total energy constraint. However, from 
Proposition]^ if v = ai, then E* has a unique solution only 
under the assumption that each AP corresponds to a different 
ai, Vi G Af. This is because E*, Vj G A/'\{i} is unique when 
V aj. If there is at least another AP, say the fcth AP, k G 
A/^\{i}, and v = ak = ai, then the optimality can be achieved 
with multiple solutions. Based on the proof of Proposition 
the solutions of Ef and E^ can be any combination, which 
is subject to that 0 < E* < £’■*’", 0 < E'^ < E’f'^ and 

TP-k I 77'-*- _ TPtot 77’-* 

^3 ■ 

Next, we propose a water-filling based distributed algorithm 
to obtain E* by finding the constant i/. Note that we only 
consider the case that each AP corresponds to a different ai, 
which can be readily modified to adapt to the case that multiple 
APs have the same ai. For convenience and without loss of 
generality, we assume that ai > a-i > • • • > a^v- As z/ is 
unique and > 0, a possible value of v could be ai,'i% G Af, 
or between an interval {ai+i,ai), where um+i = 0. Note 
that when the optimal v falls into an interval between any 
two adjacent afs, we need to apply the iterative water-filling 
algorithm | |3^ to achieve it. Specifically, after the PB releases 
the current energy price in each round, all the APs update 
their bids in parallel and feedback them to the PB within the 
current round. 

We now elaborate the steps in Algorithm [T] Each AP first 
provides its ai to the PB, then the PB constructs a descending 
sequence {ai\i^j^, sets i — 1 and announces which is 
equal to each ai, sequentially. Then, with Af the jth AP, 
Vj G A/'\{i}, computes the response, here denoted by e'^^\ 
according to (22 1 , whil^ the ith AP responds that = 


Ef 


Then, by comparing the aggregate of all the responses. 


computed by Eagg = Yi 


E^'> + EY, with El°f the 




PB’s energy will be exhausted to enhance the social welfare. 
The optimal solution E* to the problem (21 1 can be presented 
by the following proposition. 

Proposition 4. The optimal solution to the problem 0 is 
given by E* = [E*, • • • , E*ffY 


( 22 ) 


PB can decide the value of v as listed in Algorithm [1] Note 
that if V exists in an interval, then a bisection method pT] can 
be employed to find the unique value. The implementation of 
the bisection method is simple and thus omitted. It can be 
observed that Algorithm [T] can always find the exact value or 
the existing interval of v with at most N iterations, since v is 
unique and 0 < z/ < maxig 7 \/{ai}. 

According to Algorithm]^ each AP first needs to measure 
the (equivalent) channel power gains Gi and Ki, and acquire 
the transmit power pb from the PB. The values calculated by 
each AP and forwarded to the PB should be ai given in 


where u > 0 is a constant chosen to meet the total energy 
constraint Yi^Af E* = E^ and 7 i {v) is a strictly decreasing 
function on 0 < u < ai given by 


li {v) = 


Pb 


( 4 - 1 ) 


0-2 -I- Gi-q {piG^ + PbKi) 


(23) 


^According to Proposition]^ given u = oti, the response energy for the 2 th 
AP is that E* = ' However, this happens only when 

u is found to meet the total energy constraint. As proved in Appendix C, 
given u = ai, under the assumption that the 2 th AP does not know the total 
energy constraint, the possible value of E* exists in an interval [O, . 

Therefore, in this case, the 2 th AP is set to respond that E^'^^ = E^^"^ to the 
PB. Then the PB knows the interval of the 2 th AP, which will help it to make 
decisions in the following steps. 












7 



(a) Ti {Ei) versus Ei. 


(b) Si {Ei) versus Ei 


(c) VtSj {Ei) versus Ei 


Fig. 3. Graphical interpretations for the functions (n), (14} and ID in the proposed cooperative scenario. 


Algorithm 1 Water-filling based Distributed Algorithm 
1: Each AP reflects its at to the PB. 

2: The PB constructs a descending sequence {ofijigAA, sets 
i = 1 and repeats: 

1) The PB announces = Ui to all APs. 

2) The jth AP responds an optimal to the PB 
according to Proposition]^ Vj S A/^\{i}, while the 
jth AP responds 

3) The PB computes that Eagg = 
and compares Eagg with El°*: 

a) If Eagg < El°^, then set j = i + 1 and continue. 

b) Elseif E^i^gg > El°* and Ej^Jg - Ef-^ < El°\ 

then V = ai, E* = Ef\ and E* = — 

c) Else, the unique v S [ai,ai-i), which can be 
readily found via a bisection method. 


E: 


lim 


defined in (12i and the optimal response obtained by 


(22 1 . With these signals from the APs, the PB announces 


i) 


at the end of each iteration and the allocated energy E* at last 
to each AP. 

Einally, with the value of v achieved via Algorith m [T| 
and (22 1 , we can obtain E*. Subsequently, based on (111 
given in Proposition we have r* = Ti (E* ) and r* = 
K.--- . Therefore, we have already found the optimal 

solution of the original problem (|^. 


IV. Non-cooperative Scenario 

This section investigates the non-cooperative scenario, 
where the multiple APs and the PB are assumed to be 
rational and self-interested. We adopt an auction model to 
characterize the conflicting interactions among the APs and the 
PB, where the APs are bidders and the PB is the auctioneer. 
In particular, we employ a well-known Ausubel auction HD. 
which constructs an ascending-bid version of the Vickrey- 
Clarke-Groves (VCG) auction and induces truthful bidding as 
well as achieves the maximum social welfare, i.e., the global 
optimum |21 |-p4}. In the sequel, we first define the utility 
functions of both APs and the PB. Then we formulate an 
auction game, in which an auction based distributed algorithm 


“Truthful bidding means that reporting true optimal demand at every 
iteration is a mutually best response for each bidder 1221 . 


is proposed. At last, we analyze the formulated game and 
prove convergence of the proposed algorithm. 


A. Utility functions 

We first present the utility functions of APs and the PB, 
respectively. 

1) The Utility Function of Each AP: Given the unit price 
of the PB’s energy, denoted by /i, the payment from the ith 
AP can be written by 


Ei = ^Ei. 


(25) 


Therefore, with (j^ and (251, the utility function of the ith AP 
can be defined as 


(Tz, Ei^ — XiRi {Ti^ Ei'^ 

=A,(l-r,)lE 

/ Gil] {npiGi + EiKi)\ 

" V + -j 


pEi. 

(26) 


2) The Utility Function of the PB: Recall that E = 
[El, • • • ,Ejg]^. Then, the utility of the PB can be expressed 
as 

{p, E) = ^ Ti = p ^ Ei. ^27) 

zGAZ' zgW 

Note that, there is a reserve price for the PB, denoted by 
pfE) _ It is assumed that if the outcome price of an auction 
game is larger than the PB would benefit from the trade. 
Otherwise, it would not participate in the trade |24|, |25|. 


B. Auction Game 


Now, we are ready to formulate an auction game, where the 
APs are bidders and the PB is the auctioneer. Each AP submits 
bid to compete for the PB’s energy, in order to maximize its 


utility defined in (26 1 . On the other hand, the PB aims to 
maximize its revenue in ( [27] l by increasing its energy’s unit 
price. In particular, the auctioneer first announces an initial 
price p = p^^\ and the bidders respond to the auctioneer 
with their optimal demands, i.e., the bids. Then, the auctioneer 
raises the price p until the aggregate demands meet the total 
energy constraint, and meanwhile, the auctioneer concludes the 
auction and decides the final allocated energy to each bidder. 



























Algorithm 2 Auction Based Distributed Algorithm 
1: Given El°*, price step A > 0, and t = 0, the PB initializes 
the asking price 

2: With each AP computes and submits its optimal bid 


^ with (129 1 . 


by (|32|i, sets t = t + 1. 


and repeats: 

a) The PB announces to all APs. 

b) Given each AP updates its optimal bid 


obtained by (29 1 . 


c) The PB sums up all bids that 


E, 


it) 


(33 I. Then, the PB allocates E* = C. 


(T) 


the *th AP. 


and r*, the utility of the ith AP in (26i becomes to 
(Ti) = Ai(l - Ti)W 


X l0g2 1 + ^ 


GiTj [Tip^Gi + E*Ki 

(1 - Ti) 0-2 


-r*. 


We then have an optimization problem regarding Ti as 

E* 


max Ui {Ti), s.t. 0 < Ti < 1, > 


Pb 


e‘;°'> and 


The PB sums up all bids that = J^i^jv 
compares with El°^: 

1) If < El°*, then the PB ends the auction and 
quits the trade. 

2) Else, then the PB computes the cumulative clinches 


constant price ste||^ A > 0, and the initial price The 
price will be announced to all APs, and each AP will submit 
its optimal bid e\^^ based on the optimal response given by 
( [ 2 ^ . Then the PB sums up all bids that Ef'\ 

and compares with El°*. If < El°*, then the PB 
ends the auction and quits the trade. Otherwise, the PB first 
computes a cumulative clinch ED for each AP, which is the 
amount of the PB’s energy that each AP is guaranteed to be 
allocated. The cumulative clinch for the *th AP at the round 
f > 0 is given by 


cY = max ^ 0, EI°* 


E. 


- E 


it) 


(32) 


and compares with El°*: 

i) If f > El°*, then the PB records the cu¬ 
mulative clinches and sequentially 

sets + A, t = t + 1. The auction 

continues. 

ii) Else, then by setting T = t, the PB concludes 
the auction and computes {C\ }i^j\r with 


Then, the PB sets + A, f = f + 1 and updates 

to all APs. 

With the updated pG) ^ each AP submits its optimal bid 

obtained by 


E^'^ based on the optimal response ^t^*^ , Ef^ 

(29l. By comparing the aggregate bids £G) A YieJ^ 


it) 


to 


with El°\ if > El°*, then the PB records the cumulative 
clinches {cl*^}i^j\f and sequentially sets + A, 

t = t + 1. The auction continues till < El°*. By setting 
T = t, the PB concludes the auction and computes the 
cumulative clinches {C) }i^j\f according to the proportional 


In the formulated game, at each round f > 0, with a given 
price pG)^ the ith AP has an optimization problem as 

max Ui {ri.Ei^pf'*^^ , s.t. 0 < n < 1, 0 < < TiPb- 

(28) 

We denote by the optimal solution of (281. That is 

= argmax Ui (r^, Ei, . (29) 

Thus the bid of the ith AP is given We also denote 
by E* and E* the final allocated energy and the final payment 
of the *th AP when the PB concludes the auction. With E* 


rationing rule (PRR) 


cP =eP 


by 


(T-l) _ (D 




eAf 


E. 


(T-l) 


-E. 


eJ^ 


E. 


(T) 


(33) 


X Er* - ^ eP 

V i&N' / 


where cP is actually the finally allocated energy to the *th 
AP and it is a sum of its last energy bid and a proportion 
of the remaining energy (^El°* — YieAf Then, the PB 

allocates E* — cP to the ith AP. With E*, the *th AP can 
decide its final harvesting time r*. Hence, the utility of the 
ith AP can be expressed by 

u,(T*,Et) =Yii-Tnw 

GiV {T*piGi + E*Ki 


(34) 


X l0g2 1 + ' 


a-rt)<T^ 


- r* 


(30) 


where the payment from the zth AP T* is given by 

T 


n=+>cf'+Y,+>{+-cl 


(t-l) 


f=i 


(35) 


(31) Notice that the term —P 


We use T* to denote the optimal solution to problem ( |3T) , 
which is the optimal harvesting time of the ith source from 
its own AP. Analogue to the proof of Proposition [ p w e can 
readily get f* = Ti where (•) is given by ( 111 . 

Algorithm |p elaborates the auction steps. Before starting 
the auction, the PB sets up the iteration index f = 0, the 


in (351 is actually the 


amount of energy that the ith AP is guaranteed to be allocated 
with the announced price Their product will be the 

corresponding payment of the ith AP at the current round. By 

^Note that the constant step normally introduces errors to the convergence 
point |38| . However, this is not a problem for the considered auction mech¬ 
anism with the proportional rationing rule (PRR) (39| given in l |331 , which 
guarantees that the total energy of the PB can be entirely allocated at last. 
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accumulating these payments generated at each round (from 
t = 0 to t = T), we can have achieve total payment of the *th 
AP given in (35i. The utility of the PB thus is written by 

^^( 0 )^( 0 ) 


ieAf i&M 






t=l 


(36) 


C. Analysis of the Formulated Game 

In this subsection, we analyze the formulated auction game. 
We hrst derive the optimal solution to problem ( p 8 ] l given 
in the following proposition. 


Proposition 5. Given a price by the PB, the optimal 
solution to problem \28\ can be expressed as 




(^i n \ ,,(t) \ lim 


Pb 


7 , j , if /r(‘) < 

(37) 

where is defined as the maximu m a cceptable price for 
the ith AP, which is equal to at in {17), z| and 7i (•) are 
expressed in 0 and (|23p, respectively. 


Proof: See Appendix D. ■ 

Next, we prove that the proposed auction based distributed 
algorithm has the following convergence property. 


Proposition 6. The proposed auction based distributed algo¬ 
rithm converges within a finite number of iterations. 

Proof : Recall that '^he auction 

concludes when Accoring to Proposition]^ if 

the asking price > /r*®"*, then the ith AP will drop from 

the auction with = 0. If the asking price p^*'^ < /r ■®”®, then 

the ith AP always submits a non-zero bid = ji (p^*^). 
By referring to Proposition]^ we can show that 7 ^ is a 

strictly decreasing function when 0 < < p*®”® and when 

p(f) approaches p*®”®, we can get 


lim 

P,-bp.i 


(38) 


We denote by n the number of iterations when the proposed 
auction concludes. Recall that the price step is A. Thus p*^^) = 
A X n + p^^\ When n grows, p^'^'> increases. Since is 
a decreasing function of p^'^\ the summation of bids £^'^^ = 
E^'^ will decrease as p^^^ increases. Therefore, there 
always exists an that satisfies £^'^^ < El°* such that the 
PB concludes the auction. Hence, the number of iterations n 
is finite, which completes this proof. ■ 

Last, we analyze the required signaling for the computations 
in the non-cooperative scenario. According to Algorithm 
each AP first needs to measure the (equivalent) channel power 
gains Gi and Ki, and acquire the transmit power pt from the 
PB. The only value calculated by each AP and forwarded to 
the PB is the optimal bid E^*'^ obtained by (37 1 . With all the 
bids from the APs, the PB announces the price p^‘) at the end 
of each iteration and the allocated energy the allocated energy 
Et at last to each AP. 


TABLE I 

System Parameters 


I'he bandwidth W 

lOOKHz 

The noise power 

—80dBm 

The path-loss factor Q 

2 

The transmit power of the PB 

2Watt 

The transmit power of each AP pi 

IWatt 

The gain per unit throughput for each AP Xi 

lU/Mbps 

The distance between the 2 th source and its AP 

lOm 

The distance between the fth source and the FB dps^ 

10 m 

The energy conversion efficiency 77 

0.5 

The number of antennas at the FB 

4 

The reserve price for the PB 

0.001 

The price step A 

om 



(a) Convergence of the water-filling based distributed algorithm. 



-The 1st AP-Source Pair 


-The 2nd AP-Source Pair 

. 

— The 3rd AP-Source Pair , 

Approximately equal to £' 3 '" = 0.3299 
~ - -yE Auction ends and t = T 

Y 

- 

! Happens when > 

-^- 

dr 


Iteration Index t 


(b) Convergence of the auction based distributed algorithm. 
Fig. 4. The convergence properties of both the proposed algorithms. 


V. Numerical Results 


In this section, we present numerical results to illustrate and 
compare the performance of both the proposed cooperative and 
non-cooperative scenarios. We denote by dAiSi dps- the 
distance between the *th source and its AP, and the distance 
between the ith source and the PB, respectively, Vi G Jf. 
We assume that all the channels experience quasi-static flat 
Rayleigh fading and adopt a distance-dependent pass loss 
model such as = 10 “^ where dxY denotes 

dAiSi or dpSi and ^ G [2,5] is the path-loss factor. Notice 
that a 30dB average signal power attenuation is assumed at 
a reference distance of Im in the above channel model HI- 
The system parameters used in the following simulations are 
listed in Table U 

The curves in Figs. |4||^ correspond to a network setup con¬ 
sisting of three AP-source pairs with one randomly generated 
channel realization, where the channel power gains between 
the APs and sources and the equivalent channel power gains 
between the PB and the sources are 




0.0446 

0.1569 

0.8628 


X10-^ K = 


0.1616 

0.6486 

0.4379 


X 10"^ (39) 


With these parameters, we thus can compute [ai,a 2 :af\ = 
[/i'l®®®®, plf^, p^f'^] = [5.6834, 4.7802, 0.4543], [l;]®®®®, L;^®”®, 
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(a) The optimal allocated energy to each AP-source pair versus 
in the case of = 3 for both the proposed scenarios. 


(b) The optimal DL WET time to each source versus in the case 
of W = 3 for both the proposed scenarios. 


Fig. 5. The optimal allocated energy and optimal DL WET time for each AP-source pair versus E^°^ in the case of W = 3 for both the proposed cooperative 
and non-cooperative scenarios. 



- Cooperative Scenario, E[E*] 

- Non-Cooperative Scenario, E[£'*] 


(Joule) 



- Cooperative Scenario, E[rf] 

-Non-Cooperative Scenario, E[Tf] 


0.65 


0.55 


! 0.45 


i 0.35 


0.25 


3 4 5 

Ef (Joule) 


(a) The average allocated energy from the PB to each AP-source pair (b) The average DL WET time to each source versus E^°^ for both 
versus for both the proposed scenarios. the proposed scenarios. 

Fig. 6. The average allocated energy and average DL WET time for each AP-source pair versus E^^^ in the case of both the proposed cooperative and 
non-cooperative scenaiios. 


E™] = [0.1676, 0.0989, 0.3299] and = [0.6325, 

0.8307, 1.3247] based on ([T^, @ and respectively. 

Fig. g illustrates the convergence properties in the case 
= IJoul^ for both the proposed water-hlling based 
distributed algorithm in Algorithm^and the auction based dis¬ 
tributed algorithm in Algorithm g It is shown in Fig. 4(a) that 
the value of ly starts from ly = maxig 7 v/{ai} = ai = 5.6834, 
turns to ly = a 2 = 4.7802 and ly = = 0.4543. Then 

it is confirmed that ly exists between as and a 2 - With the 
bisection method, the proposed water-hlling based distributed 

®We consider T = Is to guai'antee the consistency of all physical units. 


algorithm converges w ithin few iterations and the desired 
ly is achieved. In Fig. 4(b) with the increasing price 


(as the iteration index t increases), the third AP-source pair 
hrst quits the auction when > /ig™ and its last bid 
is approximately equal to iSg*™ = 0.3299 (when the price 
step A is sufficiently small). When the price continues 
to rise, the bids of the hrst and the second AP-source pairs 
will decrease till the summation of the bids is smaller than 
Then the PB concludes the auction, sets t = T and 
allocates the energy to each AP-source pair based on ( [3^ , 
which validates the convergence of the proposed auction based 
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Fig. 7. The average social welfare (i.e., the weighted sum throughput) of both 
the proposed cooperative and non-cooperative scenarios versus in the 

case of different number of participating AP-source pairs. The social welfare 
in the non-cooperative scenario is actually the aggregate of the utilities of all 
the AP-source pairs and the PB. 


distributed algorithm. 

Fig. 1^ shows the optimal allocated energy and optimal 
DL WET time for each AP-source pair versus El°* in the 
case of iV = 3 for both the proposed cooperative and non- 
cooperative scenarios. It can be observed in Fig. 5(a) that the 
first AP-source pair, which has the largest ai = 5.6834, is first 
allocated energy by the PB. Then the second and the third AP- 
source pairs with a 2 = 4.7802 and as = 0.4543 are allocated 
energy by the PB subsequently. With the increasing of 
the optimal allocated energy E* for each AP-source pair in 
the cooperative scenario converges to E° and equals E° when 
El°^ > Meanwhile, as expected, the hnal allocated 

energy E* for each AP-source pair in the non-cooperative 
scenario can efficiently match that in the cooperative scenario 
(i.e., E*), when is small. But E* plummets to zero 
since the value of El°* exceeds that of eI°\ which 

is due to the quitting of the PB in the adopted auction 
mechanism. Note that when is sufficiently small, the 
energy payment of each AP to the PB can be ignored in its 
own utility function. In this case, the non-cooperative scenario 
can be approximated to the cooperative scenario. Thus, the 
optimal bid from each AP should approximately equal 


to E°, which leads to ~ Therefore, 

the approach of E* to E° and the plummeting of E* to zero 


occur almost simultaneously, as depicted in Fig. 5(a) 


Fig. 5(b) shows the trends of both r* and r* versus El°*. 
With the increasing of Ef^°^, the optimal DF WET time of the 
first AP-source pair first decreases and then increases while 
others remain the same. This is because the first AP-source 
pair is allocated energy first. When is small, the harvested 
energy for the source mainly comes from its associated AP, 
which corresponds to the case n > t[ (see Proposition [^. 
But, with the assistance of the PB, the AP can shorten its DF 
WET time, which leads to the decrease of t* (or r*). When 


El°* continually increases, the source can harvest more energy 
from the PB, which corresponds to the case that the DF WET 
time is dominated by the PB, i.e., = t' (see Proposition [H. 

In this case, the value of r* (or r*) grows as the value of 
increases. Moreover, we can observe the same tendency of the 
second and the third AP-source pairs when they are allocated 
energy. When > YzieAf E°, T* of each AP-source pair 
in the cooperative scenario becomes constant as E* = E°, 
almost at the same time, t* for each AP-source pair in the 
non-cooperative scenario equals that in the case each source 
only harvests energy from its own AP, which is caused by the 
quitting of the PB. 


Now we show the average performance of both the proposed 
scenarios and each curve hereafter is averaged over 10^ 
randomly generated channel realizations. Fig. demonstrates 
the average allocated energy and average DF WET time 
for each AP-source pair versus El°* in both the proposed 
cooperative and non-cooperative scenarios. E [E*], E[i?*], 
E [t*] and E [r*] are averaged for each AP-source pair and 
averaged over 10^ channel realizations. It can be observed in 
Fig. 6(a)| that the optimal allocated energy in the cooperative 
scenario E [E*] rises with the increasing of El°*. But E [E*] 
decreases with the increment of the number of participating 
AP-source pairs with a given El°* because E [E*] is averaged 
for each AP-source pair. All E [£'*]s corresponding to different 
N converge to a constant value with the increasing of El°^, 
because the larger El°^ is, the larger is the probability of 

^tot 




E° in each simulation block with one channel 
realization. For the non-cooperative scenario, E[£'*] closely 
matches E [E*] when is small, but converges to zero 
whe n El° * is sufficiently large due to the quitting of the PB. 
Fig. 6(b) depicts both E [t*] and E [r*] versus El°*. When 
El°^ increases, both E [t*] and E [r*] decrease hrst and then 
increase, and they converge to different values for the reason 
shown in Fig. 5(b) Besides, when the number of participating 
AP-source pairs rises, the turning point of the corresponding 
curve appears later. The reason is that the more AP-source 
pairs, the less energy allocated to each source from the PB. 
Then, a larger value of is needed to make the DF WET 
time at each source dominated by the PB. 


Fig. 1^ presents the average social welfare of both the 
proposed cooperative and non-cooperative scenarios versus 
El°* in the case of different number of participating AP-source 
pairs. The social welfare in the non-cooperative scenario is 
assumed to be the aggregate of the utilities of all the AP- 
source pairs and the PB. It can be observed that the network 
setup with a larger N in either proposed scenario achieves a 
better social welfare performance due to the energy resource 
brought by the new joined AP-source pairs. The average social 
welfare of the cooperative scenario increases as El°* grows 
and approaches a constant value when El°* is large enough. 
This is because when El°* > shown in Fig. 

the allocated energy to each source will be a constant value 
and independent of El°*, which leads to the saturation of the 
social welfare performance. In the non-cooperative scenario, 
as El°^ increases, the average social welfare performance 
closely matches that in the cooperative scenario when 


5(a) 
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is small, but then deteriorates due to the quitting of the PB 
when El°^ > ■ Notice that when = 0, i.e., 

without the PB, the average social welfare performance is 
always the worst. 


VI. Conclusions 

This paper investigated the joint time and energy allo¬ 
cation of a power beacon (PB)-assisted wireless-powered 
communication network (WPCN) with multiple access point 
(AP)-source pairs and a PB. We consider both cooperative 
and non-cooperative scenarios corresponding to the situations 
whether the PB provides wireless charing service to each 
AP-source pair for free or not. Moreover, the social welfare 
was maximized in the proposed cooperative scenario and the 
respective utility of each AP and the PB was maximized in 
the non-cooperative scenario. The numerical results validated 
the convergence of both the proposed water-hlling based and 
auction based distributed algorithms. It was demonstrated that 
the average social welfare of the cooperative scenario improves 
as either the number of participating AP-source pairs or the 
total energy of the PB increases, but saturates when the total 
energy of the PB is sufficiently large. Moreover, the average 
social welfare performance of the non-cooperative scenario 
closely matches that in the cooperative scenario when the 
total energy of the PB is small, but deteriorates when the total 
energy of the PB is sufficiently large, which is caused by the 
quitting of the PB in the adopted auction mechanism. 


solution to the equation (41 1 . Then, we have 1 -|—■ 
Rearranging this equality, we get 


^ 4 -1 - 


zj — 1 + Ai 


(43) 


Recall the constraints on Ti in the optimization problem (j^ 
that 0 < Ti < 1 and > —, we have 

Pb 

Ti (E,) = max jr/, , (44) 

since it is easy to check that rj < 1 and rj > 0 when Ei = 0. 
To further expand (44i, we compare the two terms 4 
Ei/pb and obtain that rj < Ei/pb only if 


Er > 


Pb 


(4-1) 


(^zj - 1^ 0-2 _|_ (^p.Q. +pi^Ki) 


(45) 


which is obtained by substituting the expression of Ai and 
Bi into ( |43] l. Therefore, the optimal solution Ti (Ei) to the 
problem can be further expressed as 


Ti (Ei) = 


4, if 0 < 

4 r (46) 

f if El^^<E,<pb, 


where is dehned as the right-hand side of (45 i. This 

completes the proof 


Appendix A 

Proof of Proposition[T] 

To proceed, we hrst derive the second-order derivative of 
Si {Ti) with respect to Ti. After some simplifications, we have 


d'^Si (n) 
dn^ 


X,W{A + B,y^ 


{1 - Ti) {1 - Ti + A^Ti + B^) ln2 


Appendix B 

Proof of Lemma[2]and Proposition[2] 


When 0 < Ei < by replacing Ti in (lOi with {Ei), 
with careful simplihcation, we have 


(40) 5, {Ei,) =A,IV 


where Ai = ^'^ 2 ' and Bi — Since ^ ^^.4^ < 0, 

Si {Ti) is a concave function of r^. Thus, the optimal solution 
can be obtained by setting the first-order derivative of Si{Ti) 
equal to zero and comparing the obtained stationary points 
with the constraints. That is, ^*^4^ = 0. After some algebraic 
manipulations, we have 


GiV {piGi -l- EiKi) 


l0g2 ( 4 ) 




(4 - 1) 0-2 -f G^ppi 

G,p {p,G 2 + E,Ki) (4 - 1 ) + Gl 




t2 J_ 

KWG^p {p.G^ + E^Ki) 


(4 - 1 ) + Gfppi 


z]a^ In 2 


(47) 


Zi \ n { zi ) — Zi + I = Ai , (41) 

where = 1 -|- Note that > 1 as > 0 and 

0 < Ti < 1. With reference to p^ , / (z^) = z^ In (z^) — z^ -f 1 
is strictly increasing when z^ > 1 and / (1) = 0. We then can 
deduce that there exists a unique solution, denoted by 4 > 1’ 
to the equality f{zi) = Ai. After some manipulations, we get 


where the second equality is based on that z| In ^z| j —zj+l = 

*^4^' (please refer to the proof in Appendix A). Thus, the 
gradient of Si {Ei) is a constant, denoted by ai, given by 


CXi — 


dS^ {E,) 
BE, 


\A¥G,pK, 

z| 0-2 In 2 


(48) 


z] = exp 



/ A,-l 

Vexp(l) 



(42) 


When < Ei < Pb, by inserting the expression of Ti 

into ([T0|l, we have 


in which W {x) is the Lambert W function that is the solution c I'c’t \ w 4 \ ^ 4 ^iEi \ 

to the equality x = Wexp(>V). Let 4 denote the optimal t) i) — i J og 2 ’ 
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where X = 


^ GiT]{piGi+pbKi) 


. By deriving the first-order 


P^ {E,) = = - -— log 2 ( 1 -f 


Pb 


Pb - E^ 


X,WX,, 


{pb - E^ + X^Ei) In 2 


Moreover, we derive the second-order derivative of Si (Ei) 
and obtain 


d^S, {Ei) 


KWX^Pb 


dE,^ {pb -E,+ X,E,f {pb -E,)ln2 


zjpb In 2 

— Oii . 


C = Y,S^{E.)-uiJ2^— 

ieX \ieX 

Taking the stationarity condition, we have 

XSi {Ei) -1^ = 0, 

and the KKT conditions are 





= 0 , 


, / , ^ b'Ph In 2 

Zi In {zi) + f ^ ^ - 1 ] z^ + I = X^, 


where AT,- = GMp.Gi+p,K,) 


= 1 -h 


derivative of Si {Ei), we have the gradient of Si {Ei) when 
j^hm < denoted by /?i {Ei), as 


XjE* 

Pb-E* 


and Zi > 1 


as Xi > 0 and Ef-"^ < E* < E°. We denote by z| > 1 
the solution of (56 1 . Then, we denote E* = ji {v) and obtain 
7 * {v) by solving z| = 1 -f as 


(50) 


li {v) = 


Pb 


(4-1) 


(z| - l) 0-2 {piG^+PbKi) 


(57) 


We then prove that zl is unique. We first define a function 
that g{zi) = Ziln(zi) -1- {Yi — 1 ) Zi -f 1 , which is the left- 


<0. (51) 


hand side of (|56| and Yi = . By deriving the first-order 


Thus, when < Ei < pb, Pi {Ei) is strictly decreasing. 

Furthermore, by inserting the expression of E)'^'^, given by 
( [T^ , into Pi {Ei), with simplification, we have 

/?* {Xb - (4 In (4) - 4 +1)) 

z • in 2 ^ ^ ^ ^ ' 

XiW f Gip {piGi + pbKP) G^ppi 


derivative of g[Zi), we have 
dg {zj) 

dzi 


= Inzi -I- Yi. 


(58) 


Therefore, as li > 0, p (zi) is monotonically increasing when 
Zi > 1. When Zi = 1, we have that g { 1 ) = Yi . Consequently, 
there exists a unique solution 4> 1 of the function g { zi ) = 
Xi, if Xi > Yi. By comparing Xi and li, we have 


(52) 


X,, -Y = 


> 


Gip {piGi + pbKi) vpb In 2 

X,W 

GippbKi vpb\n2 


X,W 


(59) 


Hence, on 0 < Fii < pb. Si {Ei) is differentiable and the 
gradient of Si {Ei) is continuous. 

Appendix C 

Proof of Proposition|4] 

According to Proposition the Karush-Kuhn-Tucker 
(KKT) conditions are both necessary and sufficient for the 


^Pblii2 /XiWG^pK^ 
X,W I 0-2 ln2 


> 0 . 


Hence, z“ > 1 is the unique solution of (56 1 . 


optimality of the problem (21 1 . To proceed, we first discuss 
the solution of a problem without the condition 0 p- E p- E°. 
By removing Q p E p E°, we have the Lagrangian as 


Furthermore, as Pi {Ei) is strictly decreasing, pi {e) is 
strictly decreasing as well. Because of the continuous property 
of XSi {Ei), we have 7 * (0) = E° and 7 ^ {ap = Ef^. 

Case 2: When Q<E^< Fif™, V5i {Ei) = a,. Thus, given 
V = Gi, any value of Ei between 0 and Ef'^ could satisfy 


(53) 

(54) 

(55a) 

(55b) 


(54i. However, by considering the KKT conditions, we have 

thatF;* = F;r‘-E,eA4\w^^;- 

The remaining case of > ai, leads to E* = 0. In 
summary, we have proved Proposition where 1 / is chosen 
to meet the total energy constraint that J^i^X Et = E^K 

Appendix D 

Proof of Proposition[5] 


The problem (28 1 can be solved by following a similar pro¬ 
cedure for problem ([^. We first apply the optimal relationship 
between Ti and Ei and simply problem ( [28] l to the following 
one 


max U., 

Ei 


Recall that the gradient of Si {Ei) is subject to that 0 < 
V Si {Ei) < ai when Q < Ei < Ef, on that basis, we begin 
to discuss the following cases. 

Case 1: When < E, < Ef, XSi {Ei) = P, {Ei), 
Pi {Ei) is strictly decreasing and 0 < Pi {Ei) < ai. Therefore, 
given 0 < 1 / < ai, there always exists a unique E* satisfying 
([^. By solving that P {E*) = v with ([Tsjl, we have 


where 


, s.t. 0 < Fi < Pb, 
Ui (F,,/rW) = 5, (F,) - 4‘^F„ 


(60) 


( 61 ) 


with Si {Ei) defined in (14i. 

Based on the gradient of Si {Ei) shown in ( [T^ , we can 
readily obtain the gradient of Ui {Ei, with respect to Ei, 
denoted by VUi {Ei), as follows. 


(56) (^i) = 


Oli - 


lim 


if 0 < Fi < F,' 
if Ff™ <E,<pb. 


(62) 
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As Pi {E]p'^') = ai and Pi {Ep is strictly decreasing (refer 
to Appendix B), WUi {Ei) is continuous on 0 < Ei < pb, 
remains constant when 0 < Ei < E^p'^ and strictly decreasing 
when Ep"^ < Ei < pb. Then, we can infer that when > 
ai, VUi (Ei) < 0, which results in that E^*'^ = 0. When 
< ai, as > 0 and according to Proposition ^ we 
have that E^*'^ is the unique solution of Pi [Ep — pS-*) = 0 
and Ep'^ < E^'^ < E° < pb- According to the analysis in 
Appendix C, we have Ef^ = 7 ^ In a summary, we 

have 



0 , if > p’p^, 

7 , (/r(‘)) , if p(*^ < /i'™, 


(63) 


where pp'^ = ai and it is defined as the maximum acceptable 
price for the ith AP. 

With the optimal E^*\ the ith AP thus can decide the 
optimal . The utility function of the jth AP with the given 
and p^*'^ can be reduced to 


(ri) =Ai(l - ri)PPlog2 

-7>e7. 


^ Gip (npiGi + E^*^Ki) 


(64) 


Based on the analysis in Appendix A, we can directly obtain 
that 


r« = 


pit) 


if 0 < Ej*^ < El, 


T'lim 


if Ef^ <EP’ <pb. 


(t) 


(65) 


Pb 


Note that when p^*^ < pp^, 7 ^ (p^*^) > pi and 

7 i {p’i'^) = Ep'^. We now can express the optimal solution 
to problem ([28]) in a compact form given by 


( 66 ) 

This completes the proof. 


n 

(^F^T)°^+Gfw7’ 


if > /r-' 


lim 
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